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Abstract. Although rare, massive stars, being the main sources 

of ionizing radiation, chemical enrichment and mechanical energy 

\^ • in the Galaxy, are the most important objects of the stellar pop- 

^^ I ulation. This review presents the many different aspects of the 

^^ ' main tool used to study these stars, i.e. spectroscopy. The first 

^ I part consists in an introduction on these objects and their phys- 

O . ical properties (mass, wind, evolution, relation with their envi- 






ronment). Next, the spectral behaviour of single massive stars is 
investigated, in the visible as well as in the X-ray domain. Finally, 
\G [ the last part of this paper deals with massive binaries, especially 

those exhibiting a colliding wind phenomenon. 

(N 
> 

CN ' 1. Astrophysics of Massive stars: Introduction 

VO ' 1.1. Main characteristics. A star is considered "massive" if it can 

ignite carbon burning in its core during the late stages of its evolution. 
\^ ! Such stars are the progenitors of black holes and neutron stars and 

O I have masses larger than lOM©. Since they go further into nucleosyn- 

thesis stages than stars like the Sun, these stars are the most important 
sources of chemical enrichment in galaxies. Such massive stars, of spec- 
O I tral type O, are blue and bright objects. Their luminosities amount to 

10^ — IO^Lq-. such objects are thus visible from far away in the Uni- 
ci ■ verse. In addition, their effective temperatures are larger than 30kK, 

meaning that the majority of their radiation is emitted in the ultra- 
^>< I violet (UV). Massive stars are therefore the main sources of ionizing 

Vj ' radiation in galaxies, and this explains why these stars are surrounded 

by bright nebulae that are H ll regions of ionized gas. 

The distribution of mass amongst stars follows a law of the form 
dN = KM~°'dM = r with N the number of stars in the mass inter- 
val [m,m + dm], K a constant and M the initial mass of stars. The 
parameter a is generally considered to be 2.35 (the Salpeter value). 
Consequently, the number of stars decreases when the mass increases. 
Massive stars are therefore very rare among the stellar population: for 
each star with a mass between 60 and 120 Mq there form 250 stars with 
masses between 1 and 2M0. This paucity also means that such stars 



Ph' 



c^ 



^Postdoctoral Researcher F.N.R.S. (Belgium). 

^BASED ON LECTURES GIVEN AT PADOVA UNIVERSITY IN DE- 
CEMBER 2005. 

1 



Y. NAZE 



Figure 1. R136a, the cluster at the core of 30 Doradus 
in the Large Magellanic Cloud (LMC), was once thought 
to be a single, supermassive star, (c) HST 



Figure 2. WR20a, the most massive system ever 
weighed, is an eclipsing binary ( Top) and a spectroscopic 
binary [Bottom): the masses of the components can thus 
be evaluated precisely. (G. Rauw, private communica- 
tion) 



are also generally distant: the nearest 0-type star, ( Oph, is situated 
between 417 and 509 light years; whereas the nearest Wolf-Rayet star 
(see below) belongs to the binary system 7^ Vel and is somewhere be- 
tween 740 and 975 light years from the Earth. 

Whilst the lowest possible mass of a massive star is generally well 
known, the question of the maximum mass of these stars is not set- 
tled yet, mainly because of the difficulties in measuring the stellar 
masses. Estimating the mass of an object can be done by modelling 
its spectrum, but these models are not always reliable, especially when 
approaching the limits of the parameter space. The use of a mass- 
luminosity relation has also been widely popular amongst astronomers 
but it can result in unrealistic large masses when the object is not 
spatially resolved. For example, R136a, at the core of the 30 Doradus 
nebula in the LMC, was once thought to be a single, 1OOO-25OOM0 
star. It was later found that R136a was actually composed of a dozen 
components (Weigelt et al. 1991, see also Fig. [T]). In fact, the only 
reliable method for deriving masses is to observe eclipsing binaries. 
Using Kepler's laws, one can then constrain all physical parameters 
of the system by observing the photometric eclipses and the spectro- 
scopic signature of the orbital motion. However, we may note that 
only young binaries, where no interactions have taken place, can lead 
to reliable, typical masses. The most massive stars detected so far by 
this method belong to WR20a and have masses of 82 and 83 Mq (Rauw 
et al. 2005, see Fig. [2]). In addition, another method has been used 
by Figer (2005) to derive an upper limit on the stellar mass: when 
observing the Arches cluster, he should have detected 20 to 30 stars 
of mass larger than I5OM0, but this was not the case (see Fig. [3]). He 
therefore concludes that there exists no star with a mass larger than 
15OM0. We emphasize that this upper limit is much larger than the 
actual largest mass observed in WR20a. 
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Figure 3. Top: The Arches cluster, observed here by 
HST, contains more than 2000 stars. Bottom : In this 
cluster, the existence of many stars with more than 150 
Mq had been predicted theoretically but none was found. 
©HST 

Figure 4. Excerpt from the digital atlas of Walborn & 
Fitzpatrick (1990). 

It can also be noted that the very first generation of massive stars 
(when there was no metaUin the Universe, i.e. population III objects), 
could have reached much larger masses (hundreds to lOOOM©). They 
could have given birth to intermediate-mass black holes, that might 
have become seeds for the supermassive black holes at the center of 
galaxies. These first stars should also have been very luminous and are 
therefore thought to be responsible for the re-ionisation of the Universe 
approximately one billion years after the Big Bang. Finally, they would 
have been the first celestial objects to build metals, sowing the whole 
Universe. However, all this is still very putative. 

1.2. Classification of Massive Stars. There exist two ways of clas- 
sifying 0-type stars: the spectral morphology or the determination of 
line ratios. 

The first criterion fits the general philosophy of the MKK classifica- 
tion scheme, where a sample of 'typical' standard spectra is constructed 
and serves as comparison when determining spectral classes. For O- 
type stars, the main tool of this kind is the atlas made by Walborn & 
Fitzpatrick (1990), which provides low-resolution spectra in the 4000- 
4700 A range (see Fig. H]). As 0-type stars are very hot, the lines useful 
for classification are indeed helium lines: the stronger the Hell lines 
compared to the He I lines, the earlier (i.e. the hotter) the star. The 
transition from O to B-type stars occurs when the HellA4542 line is 
extremely weak, barely detectable. Note that a morphological classi- 
fication is also used to classify evolved massive stars, i.e. Wolf-Rayet 
stars. 

The second criterion is more quantitative and uses the determina- 
tion of equivalent widths (EWs). The EW of a specific line represents 
the width of a rectangular line of the same area as the actual observed 
line. Conti & Frost (1977) and Mathys (1988,1989) have shown that 
the spectral type of an O star can be determined by comparing the 
EW of the HelA4471 line to that of the HellA4542 line (see Table [I]). 



Note that in astrophysics, all elements heavier than helium are called 'metals' 
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Table 1. Classification criteria for 0-type stars us- 
ing EW measurements (van der Huclit 1996): logPF' = 
log EWun - \0gEW4542; \ogW" = \ogEW4088 - 
log EW4U3; and logW" = logEWisss + log EWiese- 
Note tliat tlie recently defined spectral types 02 and 
03.5 (Walborn et al. 2002) are not presented here. 



Subtype 


logW 


Subtype 


logW 


Subclass 


03 


< -0.9 


07.5 


to 0.09 


I 


log W" > 0.3 


04 


-0.9 to -0.61 


08 


0.1 to 0.19 


III 


0.11 <logW" <0.3 


05 


-0.6 to -0.46 


08.5 


0.2 to 0.29 


V 


-0.2<logP^"<0.1 


05.5 


-0.45 to -0.31 


09 


0.3 to 0.44 


I 


logVF"'<5.15 


06 


-0.3 to -0.21 


09.5 


0.45 to 0.65 


III 


5.15 <logVF"'<5.4 


06.5 


-0.2 to -0.11 


09.7 


0.66 to 1. 


V 


log W" > 5.4 


07 


-0.1 to -0.01 











Attempts have been made to create a similar quantitative scale in other 
parts of the electromagnetic spectrum but the Conti-Mathys scale is 
still the most popular one. The luminosity class can also be deduced 
from similar EW measurements (see the same Table). 



1.3. Stellar Winds. It is well known that the Sun possesses a so- 
lar wind, which is notably responsible for generating polar auroras 
on Earth. The origin of this wind is linked to the existence of a so- 
lar corona. Outside the photosphere (~6000K), the temperature of 
the gas increases to reach lO^K in the corona. At such temperatures, 
the gas pressure is very high and it then expands naturally in the 
lower pressure surrounding regions. By this mechanism, the Sun loses 

In comparison, massive stars lose 10"'^ — 10~^MQyr~^ and they do 
not possess any corona. In fact, the presence of winds in luminous stars 
is linked to the large luminosity of these objects. In the atmosphere of 
such stars, there is a transfer of momentum from the photons of the 
stellar radiation field to the material surrounding the star. In other 
words, the winds are driven by the absorption of the stellar radiation 
and its subsequent scattering (see e.g. Lamers & Cassinelli 1999). If 
we consider an atom moving in a radial direction that absorbs a photon 
of frequency z/, its momentum becomes 

mvl = mVr + — . 
A moment later, the atom re-emits the photon at an angle a (see Fig. 
E]) and the momentum along the radial direction becomes 



mv" 



mv' 



hu' 



cos a 
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Figure 5. The stellar wind of massive stars is line- 
driven. Atoms and ions absorb and then re-emit the 
stellar radiation. 

Figure 6. Formation of a P Cygni profile, resulting 
from the superposition of an emission profile and an ab- 
sorption profile. The emission comes from light scattered 
into the line-of-sight from all regions of the wind (front, 
blueshifted wind and back, redshifted wind). The ab- 
sorption comes from light scattered away the line-of-sight 
by atoms between the star and the observer: in this part, 
the wind is coming towards us, and the absorption is thus 
blueshifted. 

The atom will mainly absorb photons at the frequency of specific 
lines. The frequency, in the rest frame, of such a line will be de- 
noted Uf). If the atmosphere were static, only radiation near the pho- 
tosphere would be absorbed. However, since there is a velocity gradi- 
ent in the wind, each region actually absorbs at a different frequency 
because of Doppler shift. If the atom is moving with a velocity Vr rel- 
ative to the star, source of photons, it will see stellar photons with a 
shifted frequency z/(l — Vr/c), so the absorbed photon has a frequency 
p = i/o(l + Vr/c) in the rest frame of the star. When the photon is sub- 
sequently re-emitted, an observer will see that an atom with a velocity 
of v'^ has emitted a photon of frequency z/' = i/o(l -|- vl/c). Therefore, 
the final velocity is 

< = Vr + ^^{l+ Vr/c) - ^(1 + V'jc)c0s{a). 

Replacing v'^. by the value derived above and assuming that v << c 
and huQ << mc^, we then find: 

//_ /^£o^ -cos(a)). 

r I rnc ^ \ J / 

Therefore, we can see that forward scattering (a=0) does not increase 
the momentum of the atom whereas backward scattering (a=180°) in- 
creases it by 2—. Since the scattering takes place randomly, we can 
calculate the average momentum gained by averaging over a solid angle 
of 47r : 

< ^mv >= ^^/o"(l - cos{a))27ism{a)da = ^. 

If the photons were coming from any direction, the net momentum 
gain would be zero. However, since they come only from the star, a 
radial acceleration appears, which results in a so-called beta velocity 
law: 

v{r) = v^{l - ^y 
where R^, is the stellar radius and (3 amounts to 0.8-1 for 0-type stars 
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Figure 7. Examples of spectra of Wolf- Rayet stars. (G. 
Rauw, private communication) 



and 1-2 for WR stars. The terminal velocities v^x, can reach 1000-3000 
kms~^ for massive stars, whereas the velocity of the solar wind is much 
lower, only 400-700 kms~^. 

The acceleration is mainly provided by the absorption and re-emission 
of UV photons in the lines of the most abundant metallic ions (CNO-f-Fe). 
As a consequence, acceleration is less effective and the mass-loss rate M 
is much smaller in a low metallicity environment. The energy gained 
by the metallic ions is subsequently shared with all other particles 
through interactions with them (Coulomb coupling), so that the whole 
wind expands. If all photons leaving the stars were absorbed or scat- 
tered, the momentum of the wind would equal that of the radiation, 
i.e. Mjnaxi^oo = L/c - this is called the single scattering hmit. However, 
note that a photon can be absorbed and re-emitted several times! 

The presence of the wind affects the emergent spectrum. The wind 
acceleration, especially in the UV, generates P Cygni line profiles (see 
Fig. E]). In addition, some emission lines (e.g. NlllAA4634,4641 or 
He IIA4686) will be formed in the wind, and not in the photosphere as 
most of the absorption lines. The spectra of Wolf-Rayet stars give an 
extreme example of this: they display exclusively broad emission lines 
formed in their thick stellar wind (see Fig. [7]). 

The radiation-driven process is an unstable one, so winds of massive 
stars are intrinsically unstable. Therefore, the wind is more patchy 
than a smooth ejection of matter. The mass-loss diagnostic lines, like 
Ha, thus arise in small, dense regions (called clumps): if the lines are 
analyzed with a spherically symmetric smooth wind model, the mass- 
loss rate will be overestimated. Generally, the filling factor of the wind, 
i.e. the ratio between the volume of the clumps and the total volume, 
is found to be 0.05-0.25. 

Since the winds are driven by radiation, the mechanical momentum 
of the wind should be related to the luminosity of the star. This rela- 
tion is often expressed as the "wind-momentum luminosity relation" : 



log Mvc 



Re 



0.5 



log Do + xlog(L*/Lc 



Q) 



where Dq is a function of the spectral type and the luminosity class 
of the star and x depends on the spectral type, luminosity class, and 
metallicity of the star (see Herrero 2005 and references therein). 
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Figure 8. In M17, a baby massive star accretes matter 
from a circumstellar disk (here seen in silhouette). © 
ESO 



1.4. Evolution. Although massive stars have been observed for cen- 
turies, astronomers still debate on how they actually form. In fact, 
these stars can ignite hydrogen core burning before they have reached 
their final mass: the radiative pressure of the light then emitted should 
prevent any further direct accretion. Three main theories have been 
proposed to overcome this problem. Massive stars could form through 
the coalescence of lower mass proto-stars, through accretion from a 
circumstellar disk (just as low-mass stars, but with a higher accretion 
rate), or through competitive accretion (the most massive clusters pro- 
ducing the most massive stars). Up to now, several evidences seem to 
favor the accretion model (e.g. Chini et al. 2004, see Fig. [8]). Note 
that the formation of massive binaries is not completely understood 
either: tidal capture and fragmentation of the unstable accretion disk 
have been proposed to explain their existence. 

Once they are formed, massive stars emit large quantities of radia- 
tion: they actually burn the hydrogen in their core at a much higher 
rate of nuclear reactions than lower mass objects. Therefore, although 
massive stars have more "fuel", they have a much shorter lifetime: a 
few million years, whereas the Sun will live ten billion years. Because 
of this short life, these stars do not generally have the time to wander 
far away from their birth place. Most of the massive stars are thus 
found in clusters, and the few isolated objects are generally considered 
as having formed in a cluster and been ejected later on, because of 
tidal interactions or following a supernova kick (if the star belonged to 
a binary system). 

If the formation of massive stars is not completely understood, this 
is also the case for the details of their subsequent evolution and death. 
The most popular evolution scenario is the so-called Conti-scenario: 
massive stars follow the path O^X^H-poor WN— >WC. WN and WC 
are two types of Wolf-Rayet (WR) stars, WN being nitrogen enriched 
and WC carbon-enriched; X is a phase that depends on the mass of the 
star, i.e. a red supergiant (RSG) phase if the star has M < 40 M© and 
a luminous blue variable (LBV) stage otherwise (for more details, see 
Maeder et al. 2005 and references therein). The wind already plays an 
important role during the main-sequence lifetime since it is responsible 
for steadily decreasing the mass of the star at a rate of about 10~^ — 
lO~^M0yr~^. However, the mass- loss further increases during the later 
stages of the evolution of massive stars: the mass-loss rate of typical 
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Figure 9. Left: rj Carinae is an LBV that under- 
went two eruptions in the 19th century, losing 1-2M0 
in the process. Right : WR124, a Wolf-Rayet star, is 
surrounded by thick ejecta from previous evolutionary 
stages. © HST 

Wolf-Rayet stars is 10~^ — lO~^M0yr~^ and gigantic mass ejections 
(whose actual trigger is not known) can even take place during the RSG 
or LBV phases (see Fig. ^. The separation between the photosphere 
and the very thick wind is thus less and less clear as the star evolves 
and material that was once in the convective core of the star, where the 
nuclear reactions happen, becomes exposed at the stellar surface. This 
explains the anomalous abundances observed in the spectrum of these 
stars: for example, WN stars show at the surface the results of hydrogen 
burning through the CNO process whereas WC stars expose layers 
that, since they were originally deeper inside the star, have been rather 
affected by helium burning. Finally, it is supposed that the massive 
stars end their life in a supernova explosion. Massive progenitors have 
even been proposed for the powerful gamma-ray bursts. 

1.5. Interactions with the surroundings. Stars are never com- 
pletely isolated, and the large energy output of the massive stars cer- 
tainly has a huge impact on their environment. These stars interact 
with their surroundings through their intense ionizing radiation, their 
powerful stellar winds, and eventually their final supernova explosion. 
The winds participate in the dissemination of the chemical elements 
built in the star, but they also transfer large amounts of mechanical 
momentum to the InterStellar Medium (ISM). We may note in this 
context that the total amount of energy released through the wind 
during the entire lifetime of a massive star is comparable to the energy 
of one supernova explosion: the wind contribution can thus not be ne- 
glected. By sweeping up the surrounding gas, the winds are therefore 
able to shape the ISM into bubbles. These bubbles have sizes of about 
ten lightyears when they are blown by single stars, but can reach thou- 
sands of lightyears if several of these stars act collectively (they are then 
called superbubbles and supergiant shells). This input of mechanical 
energy into the surroundings can be enough to induce the formation of 
new generations of stars. 

One of the first attempts to model these peculiar structures was pre- 
sented by Weaver et al. (1977, and references therein). This simple 
model has been so often used that it is generally considered as the 
'standard' one. In this scheme, a typical bubble consists of four re- 
gions, which are, starting from the star outwards, (1) a region where 
the wind freely expands, (2) a zone containing the shocked wind, (3) 
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Figure 10. Top: The structure of a typical wind-blown 
bubble in the energy-driven case (see e.g. Weaver et al. 
1977). Bottom : NGC6888, a bubble blown by a Wolf- 
Rayet star. © SDSU/MLO/Y. Chu et al. 



a region with the shocked ISM, and (4) the undisturbed ambient ISM 
(see Fig.[inD. 

The evolution of a typical bubble proceeds mainly through three 
stages (see e.g. Lamers & Cassinelli 1999). At first, the wind is not 
stopped by the ISM, and the bubble expands so fast that the radiative 
cooling does not have enough time to affect its evolution: this is called 
the 'adiabatic phase'. As the amount of swept-up matter increases, the 
shock velocity begins to decrease, and the cooling of the shell of shocked 
ISM finally becomes significant: the gas compresses into a very thin and 
dense shell. During that phase, the shocked wind becomes on the con- 
trary hotter and hotter, and has not enough time to cool. The pressure 
of this hot, shocked wind is so high that it now drives the expansion of 
the shell. This marks the beginning of the second phase, during which 
the bubble is generally called an 'energy-driven' one. Most observed 
bubbles are in this phase. Finally, the hot wind also cools, and in turn 
collapses into a dense shell. At this point, the wind impacts directly 
on the shell, and directly transfers momentum to it: the bubble is then 
called a 'momentum-conserving' one. 

During the second phase, the bubbles are detectable through (nearly) 
the whole electromagnetic spectrum. The stellar wind flows first at 
thousands of kms~^, but it nearly stops at the reverse shock that ex- 
pands with velocities of only a few tens kms~^. Its entire kinetic energy 
is then converted into thermal energy, heating up the shocked wind to 
millions of degrees. The interior of the bubble, filled with hot gas, 
will thus be detectable in X-rays. On the other hand, the thin shell of 
shocked ISM has cooled but is nevertheless typically at a temperature 
of ~ 10^ K, since it is still ionized by the stellar radiation: it will thus 
be detectable through the usual nebular emission lines (Ha, [O ill], ...). 
The interface layer between these two regions is detected through ab- 
sorptions in the UV range, but a gap between the hot interior and the 
dense shell can also reveal the presence of this layer (Chu et al. 2003). 

The Weaver et al. (1977) model attempts to describe the evolution 
of such a bubble. It relies on a few hypotheses: the ISM is initially at 
rest, and it has a uniform density no; the source of the constant wind 
is point-like and the pressure inside the bubble is uniform. To find the 
evolution of an energy-driven bubble, one then simply needs to apply 
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the conservation equations: 

Msheii = 4/3 IT R^no fJ' Conservation of mass 

— — ^^ — 2s2 = 4:71 R^ P Conservation of momentum 

at 

^ = Lw — P \^ — - Conservation of energy 

where E = 4/3 vr R^ P/{-f - 1) 

where R, Vexp and Mgheu are the radius, expansion velocity and mass 
of the thin shell of shocked ISM; hq is the ambient number density (in 
cm~^); 11 is the mean molecular weight in the Hll regioiu; P is the 
pressure inside the bubble; E is the internal energy and 7 is the adia- 
batic index of the hot gas (taken to be 5/3); and L^ = 0.5 Mv^ is the 
mechanical luminosity of the stellar wind. 

The solution to these three equations is 
Ve.p = — = 0.6/2/t 



(3850 7r)2/5^ 

The temperature and density inside the bubble can further be derived 
if the heat conduction from the shocked wind into the shocked ISM 
is considered (Weaver et al. 1977). Once the temperature distribution 
is known, the X-ray luminosity of the hot shocked wind can also be 
estimated (Chu et al. 1995). 

However, when confronted to observations, this model generally ex- 
periences several difficulties. The wind luminosity has to be decreased 
by an order of magnitude to account for the observed shell radius, shell 
expansion and ISM density. In addition, the predicted X-ray luminos- 
ity is generally not compatible with the observations. 

More specifically, Oey (1996) has performed numerical simulations 
of superbubbles, taking into account the evolution of their stellar con- 
tent, including the supernova explosions, and has compared these mod- 
els to the observational properties of several superbubbles. She found 
that these structures could be divided in two classes, one with very 
discrepant expansion velocities, and one whose velocity is more com- 
patible with the model. However, even these latter structures require 



We have used here the convention adopted by Weaver et al. (1977), although 
in the litterature, the mean molecular weight is often noted /iniH- 
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a drastic reduction of the input mechanical luminosity of the stars to 
reach a perfect agreement with the simulations. Dunne et al. (2001) 
also showed that superbubbles are generally brighter in the X-ray do- 
main than expected from their stellar content. Mass loading and/or 
off-center supernovae are often thought to be responsible for these dis- 
crepancies. 

The observed properties of single-star bubbles, like those detected 
around WR stars, have also been investigated and compared to the 
theoretical expectations. Garcia-Segura (1994) has extended the work 
of Weaver et al. (1977) to this specific case. He has taken into account 
the mass-loss evolution of the massive star prior to the WR stage, but 
his simulations agree more qualitatively than quantitatively with the 
observations. Again, the wind luminosity has to be decreased in or- 
der to match the observed radius, velocity and X-ray luminosity of the 
bubbles. 

Superbubbles and WR bubbles are rather complex objects, in which 
a lot of poorly known factors (e.g. the exact ISM density distribu- 
tion, the exact mass-loss history of the star) could influence the shape 
and evolution of the bubble. To get a realistic comparison between 
the so-called 'standard' theory and the observations, one should rather 
consider the most simple objects: bubbles blown by a single, main se- 
quence massive star interacting directly with the ISM. Such bubbles are 
called 'interstellar bubbles', but only a few of them are known. Naze 
et al. (2001b, 2002) have discovered several such structures in NllB, 
N180B, and N44, and their properties agree better with theoretical ex- 
pectations than in the case of WR bubbles and superbubbles, but the 
agreement is still far from perfect. 



2. Spectroscopy of Single Massive Stars 

Although many wavelength ranges (e.g. infrared and ultraviolet) 
provide important spectroscopic information, we will focus in this sec- 
tion on the results obtained in the visible and X-ray domains. 

2.1. Visible Domain. Apart from classifying the star, visible spectra 
are also used as input for modelling and variability studies. 

Modelling spectra constitutes a crucial step in astrophysics since 
it is the only way to derive intrinsic stellar properties like tempera- 
ture, mass- loss rate, gravity, chemical composition, rotation,... Even 
distances could be estimated thanks to modelling: once the wind- 
momentum luminosity relation (see above) is calibrated, the deter- 
mination of the radius, mass-loss rate and terminal velocity through 
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the modelling of the observed stellar spectrum leads to an estimate of 
the intrinsic luminosity of the star, hence its distance. Although mod- 
elling stellar spectra is not a trivial task, it has nevertheless greatly 
improved in recent years. The first models were indeed rather simple: 
they assumed plane-parallel and static atmospheres composed only of 
hydrogen and helium. This led to the first determinations of the tem- 
perature scales of massive stars. Very soon, a problem arose: the mass 
predicted by these models was systematically smaller than that pre- 
dicted by evolutionary models on the basis of the position of the star 
in the HR diagram (Herrero 2005 and references therein). This so- 
called mass discrepancy was mostly solved by new, improved models 
that include: 



• non-LTE (Local Thermodynamic Equilibrium) effects. Since 
massive stars possess a very intense radiation field, the radiative 
phenomena largely dominate over the collisional effects. 

• a spherical dynamic atmosphere. Plane-parallel approximation 
is valid only if the height of the atmosphere is small compared 
to the stellar radius. This is not the case for massive stars, 
since the optical depth of the wind can be significant out to 
several tens of stellar radii. In addition, the windy atmosphere 
of massive stars is expanding, and a photon emitted at one 
point can be absorbed much farther by another line thanks to 
the Doppler effect. 

• the line-blanketing. Metals, although rare, play an important 
role in the ionisation of the wind since they are very efficient at 
absorbing photons. This is especially true in the UV, where nu- 
merous lines are present: different atoms/ions can thus absorb 
radiation in the same frequency range. The inclusion of metals 
in the modelling thus results in a blocking of the UV flux, that 
leads (1) to a heating of the inner atmosphere (backwarming) 
and (2) to the cooling of the outer atmosphere (where the ion- 
isation, determined by the now reduced UV flux, is decreased). 

Such models have led to a new parameter scale for massive stars (see 
Martins et al. 2005, parameters reproduced in Table [2]). Due to the 
backwarming effect of the inner atmosphere where the absorption lines 
form, the same ionisation (i.e. helium ratios) can be reached with a 
lower effective temperature, and this leads to a reduction of the temper- 
ature scale by up to 8000K. On the contrary, for Wolf-Rayet stars, the 
emission lines are formed in the outer parts of the atmosphere, where 
the inclusion of the line blanketing results in a reduction of the ionisa- 
tion: higher effective temperatures are thus needed to fit the spectra. 
However, we may mention that the most recent models are still 1-D 
and stationary ones, and that the work continues to further improve 
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Table 2. Theoretical stellar parameters as a function of 
spectral types and luminosity classes, as determined by 
Martins et al. (2005). The effective temperature T^ff, 
here displayed in kK, is the temperature of a blackbody 
emitting the same amount of radiation as the star (there- 
fore L = iiiRlaT^j-j:); R^ is the stellar radius in R© and 
M^, the stellar mass in Mq. 
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Subtype 


Dwarf Stars (V) 


Giant Stars 


(III) 


Supei 


'giant Stars (I) 




Teff 


R, 


M, 


^e// 


R, 


M, 


'^eff 


R, 


M, 


03 


44.6 


13.8 


58.3 


42.9 


16.6 


58.6 


42.6 


18.5 


66.9 


04 


43.4 


12.3 


46.2 


41.5 


15.8 


48.8 


40.7 


18.9 


58.0 


05 


41.5 


11.1 


37.3 


39.5 


15.3 


41.5 


38.5 


19.5 


50.9 


05.5 


40.1 


10.6 


34.2 


38.0 


15.1 


38.9 


37.1 


19.9 


48.3 


06 


38.2 


10.2 


31.7 


36.7 


15.0 


36.4 


35.7 


20.3 


45.8 


06.5 


36.8 


9.8 


29.0 


35.6 


14.7 


33.7 


34.7 


20.7 


43.1 


07 


35.5 


9.4 


26.5 


34.6 


14.5 


31.2 


33.3 


21.1 


40.9 


07.5 


34.4 


8.9 


24.2 


33.5 


14.3 


29.1 


31.9 


21.7 


39.2 


08 


33.4 


8.5 


22.0 


32.6 


14.1 


26.9 


31.0 


22.0 


36.8 


08.5 


32.5 


8.1 


19.8 


31.7 


13.9 


24.8 


30.5 


22.2 


33.9 


09 


31.5 


7.7 


18.0 


30.7 


13.7 


23.1 


29.6 


22.6 


32.0 


09.5 


30.5 


7.4 


16.5 


30.2 


13.4 


21.0 


28.4 


23.1 


30.4 



the models. 



Another possible analysis of the spectrum of a single star is to in- 
vestigate its variability. A very efficient tool for this task is the Tem- 
poral Variance Spectrum (TVS) that was defined by FuUerton et al. 
(1996). Consider a dataset of A^ normalized spectra with the same 
wavelength sampling and arrange them in a matrix S where Sij is the 
jth wavelength element of the ith spectrum. To search for variability, 
the spectra have to be compared to a mean spectrum, and the differ- 
ences tested statistically. However, the signal-to-noise ratio is different 
for each spectrum and a weighting is therefore needed in order to per- 
form meaningful statistical tests. To this aim, we define the weighting 
factor Wi as a'^/af^, where aic is the inverse of the signal-to-noise in the 

continuum of the ith spectrum, and cTq is equal to ( -^ ^j^^ crr^ j , 

i.e. ao is the inverse of the rms signal-to-noise in the continuum of 
the dataset. The weighting factor therefore enables to reduce the im- 
portance of low quality (i.e. low signal-to-noise) spectra. It is also 
important to take into account the wavelength-to-wavelength varia- 
tions of the noise: indeed, the signal-to-noise ratio is higher in emission 
lines, but the intrinsic variations are also higher and larger deviations 
are thus normally expected for these lines. The observed deviations 
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Figure 11. Mean visible spectrum of HD 191612 (top) 
and Time Variance Spectrum (TVS, bottom). (From 
Naze et al. 2005) 

Figure 12. A corotating interaction region, with a spi- 
ral form (dashed line), can modify the spectrum of mas- 
sive stars. A component will first appear in the line of 
sight at intermediate velocities (case A, solid line). As 
the star revolves, only parts of the spiral can be seen 
at low and high velocities (case B, solid lines) and the 
observed absorption component therefore seems to mi- 
grate towards long and short wavelengths, the relative 
importance of the blueshifted vs. redshifted component 
depending on the physical parameters of the region and 
the star (e.g. rotation rate, change in the mass-loss rate 
and/or rotation). 

Figure 13. Schematic view of the 6^ OriC system. 

Figure 14. Variations of the EW of the Ha line (top), 
of the He I A 4471 and He ll A 4542 lines (middle) and Hip- 
parcos photometry of HD 191612 (bottom). (From Naze 
et al. 2005) 

have thus to be "normalized" by a correcting factor reflecting the ex- 
pected noise. If the exposure time is sufficient, the instrumental noise 
is negligible compared to the photon noise and a good correction factor 
crfjio-fc would be simply S'jj, the spectrum itself: aij is then < 1 



a 



V — ^ij 

if there is an absorption line (since a lower signal implies a lower noise) 
and > 1 in the case of an emission line above the continuum. Taking 
these correcting factors into account, the TVS at a given wavelength 
is then defined as: 

''J 

with the mean spectrum Sj = -^ X]j=i '^i^ij 
This TVS follows a distribution (TqXn-i (x^ being here the reduced chi- 
squared), and a statistical test can therefore be performed easily: the 
deviations are generally considered significant if they exceed the 99% 
level. An example of such a TVS is shown in Fig. [TTl 

Once a variability is detected and if the spectra are sufficiently nu- 
merous, the dataset can be searched for periodicities in the observed 
variations. This can be done for example with the Generalized Fourier 
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Transform (see Heck et al. 1985 and remarks in Gosset et al. 2001) 
that extends the Fourier Transform to the case of non- regular samphng. 

For massive stars, a stochastic variabihty is generally expected be- 
cause of the presence of the unstable stellar wind that generates short- 
lived small-scale structures. In fact, the instability of the wind is inti- 
mately linked to the line-driven mechanism. Consider an atom or ion at 
a distance r of the star and moving at a velocity Vr- It absorbs photons 
at a frequency z/q in its reference frame, i.e. a frequency vq{1 + Vr/c) in 
the reference frame of the star (cf. above). The velocity can be slightly 
perturbed and become v + 5v. If 5f > 0, the atom/ion will absorb 
photons of higher frequency, of which plenty are available, and it will 
therefore accelerate even more: the perturbation is thus amplified. On 
the contrary, ii 5v < 0, the atom/ion can only absorb photons of lower 
frequency but these have already been absorbed by the slower material 
closer to the star and are no more available: the particle will therefore 
decelerate even more with respect to the unperturbed Vr velocity law, 
and the perturbation is again amplified. Hence, any slight perturbation 
of the velocity is doomed to be amplified, provoking the formation of 
small-scale structures (so-called clumps) in the wind. These structures 
are thought to produce a stochastic variability of the spectrum. 

Sometimes, however, variations appear to be regular, with periods 
ranging from a few hours to several years. The most common sources 
of variability include: 

• Pulsations. It is well known that the spectra of Cepheids change 
as these stars pulsate radially. Non-radial pulsations of lower in- 
tensity can also modify the line profiles and magnitudes of stars. 
Generally, asteroseismology has focused on low-mass stars, but 
a few massive objects (e.g. ( Oph, HD 152219, HD 93521) ap- 
parently also display pulsations. 

• Structures in the Wind. Structures in the wind on a rather large 
scale can appear and disappear, modifying the observed spec- 
trum. For example, if the stellar surface harbors a cold (resp. 
hot) point, the mass-loss rate above that point will be decreased 
(resp. increased) and the velocity of the gas will become larger 
(resp. smaller) because of the reduced (resp. increased) absorp- 
tion: this modified wind will soon collide with the "normal" 
surrounding gas. Due to the stellar rotation, spiral structures 
might then appear, and give rise e.g. to Discrete Absorption 
Components (DACs, see Fig. [12] and Cranmer & Owocki 1996). 

• Magnetic Fields. Because of the very broad lines in the spec- 
tra of the massive stars, it is very difficult to estimate their 
magnetic field and so far there are only two cases, the stars 
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9^ OriC and HD 191612, where it has been measured with cer- 
tainty. In the 9^ Ori C system, the wind material is funneled 
by the magnetic field towards the magnetic equator, creating 
a dense equatorial region where some emission lines can arise. 
A recurrent modulation of the spectrum then appears because 
the magnetic axis is not aligned with the rotation axis (hence 
the name magnetic oblique rotator used for 6'^ OriC): different 
parts of the "disk" are therefore seen at different phases of the 
rotation cycle (see Fig. [13] and Stahl et al. 1996). 

In this context, the peculiar variability of the Of?p stars needs to be 
mentioned. The Of ?p category was introduced by Nolan Walborn in 
1972 to describe two stars, HD 108 and HD 148937, with spectra that 
were slightly different from those of normal Of supergiants. Notably, 
they present C ill lines around 4650A with an intensity comparable to 
that of the neighbouring Nlll lines. In addition, their spectra show 
sharp emission lines and some P Cygni profiles. A third star was soon 
added to this new class, HD 191612. The observation of HD 108, the 
best studied member of this class, led to conflicting results in the past, 
with explanations for the radial velocity variations ranging from bi- 
nary motion to stochastic wind instabilities. Using a 15yr monitoring 
campaign of the star. Naze et al. (2001a) discovered that the star 
actually underwent long-term line profile variations: the Balmer lines 
and the He I lines passed from emission or P Cygni profiles to absorp- 
tions while other emission and absorption lines were unchanged, like 
He II A 4542. These variations appear recurrent with a timescale of ap- 
proximately 50-60 years. A few years later, Walborn et al. (2003, 2004) 
reported a very similar phenomenon in the spectrum of another Of?p 
star, HD 191612 (see Fig. [T^ . but the timescale appears much shorter, 
about ~540 days. The same timescale was subsequently detected in 
Hipparcos photometry (see Naze et al. 2005 and references therein). 
Investigations to determine the exact nature of these peculiar stars are 
still ongoing. 

2.2. X-rays. Because of technical difficulties. X-ray astronomy was 
born quite recently. Indeed, rockets or satellites are needed to over- 
come the absorbing effect of our atmosphere - doing X-ray astronomy 
was thus not possible before World War II. In addition, the efficiency 
of X-ray detectors and telescopes has improved very slowly, and that 
explains why the first generation of "great observatories" (i.e. Chan- 
dra and 'XM.yi- Newton) was only launched in 1999. 

Spectroscopy in X-rays can be performed in 3 different ways: 

• CCDs. CCDs for X-rays can provide much more than a simple 
image. Due to the low luminosities of astronomical objects in 
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Figure 15. X-ray spectrum of 9Sgr, obtained with 
X.M.M.-Newton using a CCD (top) or a grating (bot- 
tom). In the low resolution spectrum, the lines seen at 
high resolution are blended, forming a bell-shape pseudo- 
continuum, (from Rauw et al. 2002) 



the X-ray range, the X-ray photons can in fact be recorded one 
at a time. The electron shower generated at the arrival of the 
photon will therefore be recorded precisely in position and also 
in intensity. As the number of electrons is directly proportional 
to the energy of the incident photon, CCDs provide a cheap 
and simple way to do spectroscopy, although only at a rather 
low-resolution {R = dE/E ~ 10 - 50). 

• Gratings. Gratings (either in transmission or reflection) can be 
used in the X-rays with only little modification compared to the 
visible domain. Such instruments provide a higher resolution 
than CCDs: /? ~ 200 - 2000. 

• Bolometers. As for CCDs, X-rays are detected once at a time in 
bolometers, where the photon energy is converted into thermal 
energy of the electrons. Since a higher frequency photon will 
lead to a larger temperature increase of the bolometer, spec- 
troscopy is a direct by-product of the use of bolometers. They 
provide high-resolution spectra {R ~ 1000), but they need to 
be cooled to O.IK (because an X-ray photon will provoke a AT 
of only a few O.OOIK!). The Japanese observatories Astro-E 
and Astro-E2 should have used the first bolometers for X-ray 
astronomy, but the former exploded after launch and the lat- 
ter lost all its liquid helium (necessary to cool the bolometer) 
shortly after orbit insertion. 

An example of spectra obtained with the first and second methods is 
shown in Fig. [151 

As far as massive stars are concerned. X-ray astronomy really began 
25 years ago. At that time, the Einstein observatory had just been 
launched and NASA was trying to calibrate it by observing well-known 
sources. The observation of one of these sources, Cyg X-3, revealed four 
nearby spots (see Fig. [16]). At first thought to be due to an instrumen- 
tal effect, the spots were soon found to mark the discovery of X-ray 
emission from 4 massive stars belonging to the Cyg 0B2 star cluster. 
Indeed, Einstein and its followers showed that X-ray emission is very 
common among massive stars. 
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Figure 16. The X-ray emission from massive stars 
was discovered serendipitously in December 1978, when 
Einstein observed the bright Cyg X-3 for cahbration 
purposes. The four "spots" above Cyg X-3 correspond 
to the massive stars Cyg 0B2 #5, 8, 9 and 12. (c) 
Einstein 



However, the exact origin of that emission is still under debate. Some 
authors had first proposed that a corona at the base of the wind, anal- 
ogous to what exists in low-mass stars, could be responsible for the 
high-energy emission of massive stars. However, several observational 
objections against such models were raised (see Owocki & Cohen 1999 
and references therein) : absence of a strong attenuation by the stellar 
wind (this suggests that the source of the X-ray emission lies signifi- 
cantly above the photosphere, at several stellar radii), too low an X-ray 
output, inconsistencies between UV and X-ray predictions compared to 
observations,... As an alternative, a scenario based on the instability 
of the line-driven mechanism has been proposed. In fact, an unstable 
line-driven wind does not flow at the same velocity everywhere and 
shocks between the different parts are expected, causing the formation 
of dense shells which will be distributed throughout the whole wind. 
At first, the forward shocks between a fast wind and the ambient slow 
( "shadowed" ) material were considered as the probable cause of X-ray 
emission, but subsequent hydrodynamical simulations rather showed 
the presence of strong reverse shocks which decelerate the fast, low 
density material. However, the resulting X-ray emission from such ma- 
terial after it crossed the reverse shock is very low, and can probably 
not explain the level of X-ray emission observed amongst O stars. More 
recent simulations by Feldmeier et al. (1997) have shown that mutual 
collisions of dense shells of gas compressed in the shocks would lead to 
substantial X-ray luminosities, comparable to the observed ones. Such 
models also predict significant short-term variations of the X-ray flux 
but since this is not observed, it was concluded that the winds are most 
probably fragmented (or clumpy), so that individual X-ray fluctuations 
are smoothed out over the whole emitting volume, leading to a rather 
constant X-ray output (Feldmeier et al. 1997). In addition, a supple- 
mentary X-ray emission may result from other mechanisms, which are 
not necessarily present in every massive star. For example, an accreting 
compact object will generally emit a wealth of X-rays and can have a 
drastic impact on the stellar wind structure (e.g. Kallman & Mc Cray 
1982). In binary systems containing two hot stars, a colliding-wind 
phenomenon (see next section) and/or inverse Compton scattering by 
relativistic electrons accelerated by the shocks (Chen & White 1991) 
will also lead to additional X-ray emission. Finally, the wind from both 
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hemispheres, deviated by a magnetic field, can colhde in the equato- 
rial regions, and provide another substantial source of X-ray emission 
(Babel & Montmerle 1997). 

On the observational side, it was soon found that the X-ray luminos- 
ity scales with the bolometric luminosity of massive stars. Although 
quite dispersed in the past, the data now show the relation to be rather 
tight: log(Lxm0.5 - 10. keV) = \og{LBOL) - 6.91 ± 0.15 (Sana et al. 
2006b). On the other hand, low-resolution spectroscopy has unveiled 
the fundamental characteristics of the X-ray emission. First, it is not of 
the blackbody-type: as the heated gas is optically transparent at these 
wavelengths, the observed emission corresponds to the superposition of 
discrete emission lines. This emission can thus be fitted by optically 
thin thermal plasma "mekal" or "Raymond-Smith" models. Moreover, 
the temperature of the emitting gas is found to be about O.SkeV, or 
6MK, and the wind absorption is generally low, except for Wolf-Rayet 
stars. 

The high-resolution spectra enable to go further into the analysis of 
the X-ray emission. Two main types of studies can be undertaken. The 
first one is related to the fir triplets. These lines are seen in helium-like 
ions and correspond to transitions from the n = 2 level to the n = 1 
ground level. The / line, or forbidden line, arises from the transition Is 
2s (^5*1) —>■ Is^ CSq) whereas the r (for resonance) line is linked to the 
transition Is 2p (^Pi) — > Is^ (^•S'o) and the i (for intercombination) line 
to the transition Is 2p (^Pi) -^ Is^ (^'S'o). In stellar coronae, the ratio 
of the / and i lines mainly depends of the electron density. However, in 
the hot plasma surrounding massive stars, the UV radiation plays an 
important role by coupling the upper levels of these two lines thereby 
reducing the / line in favor of the i line. Therefore, the 9fJ = f/i ratio 

is a diagnosis of the dilution factor 0.5 I 1 — Wl — (-^) J of the UV 

radiation, i.e. a diagnosis of the distance r from the stellar surface 
where the X-rays are emitted (see Fig. [T7I) . For example, the X-rays 
from ( Pup and 9 Sgr apparently form at a few stellar radii, as predicted 
by the standard modeo, whereas the f/i ratio observed for 9^ Ori C 
suggests that the X-ray emission arises very close to the photosphere 
(Rauw 2005b). Note that the G = ^^ ratio is sensitive to the temper- 
ature of the gas. 



If the X-ray emitting material is distributed throughout the wind, the hues 
should form near the radius where t = 1 at the specific energy of the line. This is 
the case for 9 Sgr (see right part of Fig. [T7]). 
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Figure 17. Analysis of the fir lines of 9 Sgr. Top: The- 
oretical f /i ratio for different helium-like ions with the 
90% confidence interval of the observed ratio for Neix 
(thick solid line). Bottom: Radius of optical depth unity 
as a function of wavelength: the observed ratio is com- 
patible with the radius of unity optical depth, (from 
Rauw et al. 2002) 

Figure 18. X-rays lines at high resolution. Top: ( Pup 
(from Cassinelli et al. 2001) Bottom: ( Ori (from Wal- 
dron & Cassinelh 2001) 

In addition, high-resolution spectroscopy also enables to investigate 
the detailed morphology of the X-ray lines, leading to additional physi- 
cal information. If the wind were optically thin without any absorption, 
the lines would appear flat-topped (see right part of Fig. [TSj) . However, 
the wind absorbs part of the X-rays that it has emitted itself. This ab- 
sorption will be larger if there is more material between the emission 
region and the observer, as is the case for X-rays emitted in the reced- 
ing part of the wind. The absorption will thus affect more the red part 
of the line: the observed line will therefore appear blueshifted (see right 
part of Fig. [T8l) . In addition, the wind expands with a large velocity, 
so that the lines should appear rather broad. Such broad, asymmetric 
lines are observed for ( Pup and 9 Sgr. However, ( Ori displays broad 
but symmetric lines (Fig. [THl Rauw 2005b). To explain this differ- 
ence, Oskinova et al. (2005) have proposed to consider a porous wind 
consisting of optically-thick clumps: although each clump can absorb 
efficiently the radiation, the X-rays can still escape freely by passing 
between them. Therefore, the absorption does not depend on the opac- 
ity of the wind (clumps), but on the spatial distribution of the clumps. 
If the clumps are tightly packed, the wind is nearly homogeneous and 
the lines will be skewed. On the other hand, if the clumps are rare 
and distant, the lines will be symmetric. Finally, 6 Ori presents nar- 
row symmetric lines that could be due to a colliding wind phenomenon 
(see below) while the lines of 6^ Ori C can be easily explained by the 
confined wind model (Rauw 2005b). 

3. Spectroscopy of Massive Binaries 

3.1. Colliding Winds : Introduction. As we have seen before, mas- 
sive stars blow dense and powerful stellar winds. If such stars belong 
to a binary system, a collision between the two winds is unavoidable. 
Since the winds are supersonic, the shock between them is a strong one 
and the gas will become very hot and dense after the collision. This 
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Figure 19. Hydrodynamical simulations for HD 152248 
(left, radiative case) and HD 93403 (right, adiabatic 
case). (H. Sana, private communication) 



phenomenon was predicted a few decades ago, but it was only con- 
sidered seriously since recent years, when the observational evidence 
began to accumulate. Some theoretical considerations will first be pre- 
sented, before reviewing the observational data. 

The temperature of the gas after the shock can be evaluated by T = 
-j^ (with fj, the mean weight of wind particles and v± the component 
of the pre-shock wind velocity perpendicular to the shock surface). For 
solar abundances, this equation can be written as T(K) ~ 1.410^f^g 
or A;T(keV) ~ 1.2t>^g, if v±^8 is the velocity v± expressed in thousands 
of kms~^. For typical winds {vj_ ~ 2000km s~^), the temperature can 
reach 60MK. To understand the physical properties of the post-shock 
gas, we can use the ratio between the characteristic timescale of ra- 
diative cooling and the time to escape the shock zone (Stevens et al. 

4 

1992): X = 't^ ~ ~^ — where xy is the separation between the con- 
sidered star and the stagnation point (the intersection of the contact 
surface with the axis joining the two stars' centers), expressed in lO^km 
and the mass-loss rate M should be given in lO~^M0yr~^. For values 
of X close to or larger than one, the cooling does not have the time 
to play a role and the shock can be considered as adiabatic: the gas 
remains at a high temperature. In this case, no optical emission line 
should form in the collision zone. On the contrary, when x << 1, the 
cooling is very efficient and the shock radiates a lot (i.e. emission lines 
over a broad range of wavelength, including the optical domain, will 
now be generated in the shock zone). In this case, the shock zone is 
compressed and subject to many instabilities (see Fig. [19] and Stevens 
et al. 1992). 

The geometry of the collision zone can also be derived rather easily 
since the contact surface of the two winds corresponds to the equilib- 
rium between the two wind ram pressures (Stevens et al. 1992 and Fig. 

ms: 

Pi Wi,l = P2 v\<^ 

or pi v\ COS^ 01 = P2 v\ COS^ 02 

Since the continuity equation implies M = 4 vr r^ t; p, the above relation 
can be re-written as: 

Mii^ COS^ 01 = Mli| C0S2 02 
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Figure 20. Geometry of the colliding wind region. (G. 
Rauw, private communication) 

Figure 21. Contact surfaces for different momentum 
ratios. (From Sana & Rauw 2003) 

As we have seen before, the velocity in the wind follows a /5-law, i.e. 
v{r) = foo(l — — )^ and we can therefore write: 

cos 62 = TZ— , ~„ / ,^^2 cos 61 
if we define the on-axis momentum ratio as 7^ — ^ ^'^'^ 



N 1/2 

^M2Vao,2 

The equation is indeed equivalent to: 

cos 02 = A COS (j)i y- 

where A = 7^ (1 - Ri/nf^/^/ll - R^lr^f''!'^. 

In addition, we know that 7r/2 — (pi = f3 — 61 and tt/2 — 02 = ^2 ~ /? 
(see Fig. [201), ^^^ therefore 

cos 01 = sin P cos 6'i — cos P sin 6*1 
cos 02 = cos j3 sin 6*2 — sin /3 cos 62 

The equation then becomes 

ri \ tan /3 cos 0i —sin 61 

r2 sin92— tan/3 COS62 

And finally, the equation of the contact surface is 

dz L. n (Ar|+rJ)z 

dx ^ Atj x+r'f (x—d) 

with ri = y/x^ + z^, r2 = \/{x — dY + z^ 

If the winds have reached the terminal velocity before they collide, 
A simplifies to TZ. In this specific case, the on-axis momentum ratio 
TZ can be easily physically interpreted. First, it enables to find the 
stagnation point. In fact, in this case, the two angles 0i and 02 are 
zero and IZ directly gives the ratio of the distances between the stag- 
nation point and the stars. Therefore, the stagnation is closer to the 
star with the weaker wind (see Fig. [211) . In addition, TZ also gives the 
form of the shock because the half opening angle of the shock cone (in 

°) was found empirically to be ~ 120 (1 — ^-^— ) TZ~^'^. The shock 

thus wraps around the star with the weaker wind (see Fig. [?T|) . 

In the above discussion, several effects have been neglected. For ex- 
ample, in some tight binary systems, there will be a deflection of the 
shock zone because of the orbital motion. Moreover, if the system is 
eccentric, the colliding wind (CW) phenomenon could appear only near 
periastron, i.e. when the stars are closer to each other. In addition. 
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Figure 22. Definition of tlie axes for tomograpliic anal- 
ysis (see text). 

since tlie winds are driven by tlie scattering of UV photons, the wind 
of one star is affected by the presence of the radiation field of the other 
star. This effect is called 'radiative inhibition': the radiative pressure 
from a companion being able to slow the radiation-driven winds, a 
weaker CW shock will result. Finally, the asymmetry or dumpiness 
of the winds has also been neglected and only starts to be taken into 
account in some hydrodynamical simulations. 

3.2. Signature in the Visible Range. A binary can be easily de- 
tected with spectroscopy, except if the system is seen face-on, since the 
lines regularly shift from the blue to the red side of the spectrum, in 
harmony with the orbital motion. If the moving lines of only one star 
are seen, one talks about an SBl (spectroscopic binary with one com- 
ponent detectable); if the lines of both stars are observed, the system 
is called a SB2. Normally, all lines detected in the composite spectrum 
of the system belong to one star or the other but, in the case of a ra- 
diative CW, additional emission lines appear. Since these lines are not 
formed in or close to the photosphere of the stars, they do not follow 
the orbital motion of the system. 

Doppler tomography can help to better determine the properties of 
these peculiar emissions but to apply that technique, a good spectro- 
scopic coverage of the orbital cycle is crucial. Two versions of the 
Doppler tomography are available: the simple S-wave analysis and the 
more sophisticated Doppler mapping (Rauw 2005a). Both request the 
definition of specific axes {x, y): x is the binary axis, from the primary 
to the secondary star; y is perpendicular to x, and in the direction of 
the motion of the primary star (see Fig. [22l) . These axes are not fixed 
but rather rotate with the orbital motion. If an emission component 
(either a discrete one or the position of the peak of a broad emission 
line) has fixed velocities in this reference frame, it will appear with 
a velocity v{(f)) = V{(t>) sini = — \4 cos(27r0) sini -|- V^sin(27r0) sini or 
—VxCos{27i(j)) + VySm{2n(f)) throughout the orbital cycle (0 being the 
orbital phase). The above combination of a cosine and a sine generates 
a radial velocity curve that displays an S-shape in the dynamical spec- 
trum, hence the name 'S-wave analysis'. In practice, one thus tries to 
fit an S-shape function to the observed radial velocities of the chosen 
emission component. Then a Doppler map reporting the fitted {vx,Vy) 
is created. This map generally also shows the star's position in the 
velocity space, i.e. (0, — i^i) for the primary and {0,K2) for the sec- 
ondary, considering that the system is not eccentric (the stars never 
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get closer to or farther away from each other). The velocity equivalent 
of the Roche lobes can also be displayed, as is the case in Fig. | 



To analyze broad emissions, a tool first developed for medical imag- 
ing must be used: the Radon transform. It is defined as g{s,9) = 
J_ f{j, k)dt where (s, t) and (j, k) are two orthogonal reference frames 
with an angle 6 between j and s, i.e. j = scos{6) — tsm{6) and 
k = ssin(6') + tcos(6'). In astronomical spectroscopy, g{s,9) is actually 
I{v, (f)), i.e. the intensity of a spectrum at a given velocity and a given 
phase, which results from the integration of the emissivity function 
along the line-of-sight through the {vx,Vy) space. As a consequence, 
one does not use the Radon transform to get Doppler maps, but the 
inverse Radon transform. Practically, the spectra should first be fil- 
tered to suppress the high-frequency noise which would degrade the 
Point Spread Function of the resulting map. Next, for each observed 
phase, one must find the intensity in the spectrum at a velocity corre- 
sponding to each pair {vx,Vy). Finally, the resulting intensities along 
the orbital cycle are added, using a weighting to take into account the 
different phase intervals covered by each observation (Rauw 2005a). 
If the emission corresponds preferentially to a specific region in the 
{vx, Vy) space, a peak will appear at that position of the Doppler map; 
on the contrary, the signal at other velocities appears with random in- 
tensities and will therefore cancel out (see Fig. 



While Doppler tomography can be very useful, it must be reminded 
that the Doppler maps display the velocity field. They are NOT 'usual', 
spatial maps, and should thus not be interpreted as such: components 
close in the velocity space can actually be very distant in position. An- 
other example is the Doppler map of an accretion disk: the higher ve- 
locities are reached closer to the star, therefore, the inner (resp. outer) 
part of the disk will appear on the outside (resp. inside) in the veloc- 
ity space! In addition, the tomography should be applied with care to 
eclipsing binary systems (the relative importance of the components 
would then be biased). The same holds for the systems where emis- 
sion arises outside the orbital plane (which is unfortunately the case 
for winds of massive stars). 

On the other hand, other signatures of the CW phenomenon can also 
be detected. For example, in HD 152248 (07.5III+07III), Sana et al. 
(2001) discovered that the strength and width of some emission lines 
was varying (see Fig. [221). I^ f^ct, these emission lines were broader 
and stronger at quadrature (i.e. when the stars have maximum ra- 
dial velocities) than at conjunction. This can be easily explained by 
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Figure 23. Doppler map showing the resuhs of the 
S-wave analysis for HD 149404 (07.5I+ON9.7I). The 
crosses indicate the velocity of the center of mass of the 
binary components while the dashed line is the equiva- 
lent of the Roche lobe in the velocity space. The dif- 
ferent symbols stand for different emission lines: filled 
circles = Ha; open circles = H/?; filled triangles = 
NlllAA5932,5942; open triangles = NlllAA4634-41; filled 
squares = S lvAA4486,4504; open squares = ClllA5696 
and stars = HellA4686. (From Rauw et al. 2001) 



Figure 24. Doppler tomography of the HellA4686 line 
of LSS3074 (04f"'"+06). (G. Rauw, private communica- 
tion) 



Figure 25. Full widths of the base of the emission com- 
ponent of the HellA4686 (left) and Ha (right) lines of 
HD 152248. Different symbols refer to different instru- 
ments. Note that the absorption components of the stars 
have been subtracted from the original line before mea- 
suring its width. (From Sana et al. 2001) 

the presence of a planao CW region just between the two stars. The 
emission components produced in this dense CW region would be oc- 
culted at conjunction phases: the equivalent width of the lines should 
be smaller then. In addition, since the shock is almost perpendicular 
to the axis of the system, the radial velocities of the particles escaping 
from the wind interaction region show a broader distribution when our 
line of sight is aligned with the interaction zone (i.e. at quadrature) 
than at conjunction. 



3.3. Signature in the X-rays. In view of the high post-shock tem- 
perature. X-ray emission appears as an obvious signature of CW phe- 
nomena. Such signatures have indeed been found in many binaries, 
and they are generally characterized by: 

• Large X-ray luminosity. As CW represent an additional phe- 
nomena, the X-ray luminosity of binary systems with CW should 
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since both stars have very similar stehar and wind properties 
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Figure 26. Lx - Lbol relation for the NGC6231 clus- 
ter. The two colliding wind binaries (circled) clearly lie 
above the canonical relation. Note that the X-ray emis- 
sion of HD 326329 is probably contaminated. (From Sana 
et al. 2006a) 



exceed the simple combination of the two individual X-ray lu- 
minosities of the massive stars. Therefore, when comparing X- 
ray luminosities to bolometric luminosities, CW systems appear 
above the 'classical' Lx — Lbol relation (see Fig. [26]) . 

• High temperature. The X-ray emission from massive stars gen- 
erally displays a rather low temperature kT (about O.SkeV). 
However, the emission from CW arises in hotter plasma and 
should therefore present higher kT . Non-thermal emission com- 
ing from inverse Compton scattering (see below) could also be 
observed in colliding wind binaries (CWB). 

• Modulation of the X-ray flux. As the binary system rotates, dif- 
ferent regions come into view/the line-of-sight. A modulation of 
the X-ray flux is therefore expected due to changing absorption 
of the CW emission. In eccentric systems, the variation of the 
separation d between the stars induce variations in the emitted 
X-ray flux {Lx oc v~^^'^d"^ , see Stevens et al. 1992). 

For example, in the eccentric binary system HD 93403 (05.5I-I-07V), 
modulations of the X-ray flux are observed in different energy ranges 
(see Fig. [271 and Rauw et al. 2002). The soft X-ray emission most 
likely arises in the outer regions of the individual stellar winds and the 
variability in this energy range is probably associated with opacity ef- 
fects. In the medium energy band, these effects are much smaller and 
the observed variation is consistent with a l/rf modulation, where d is 
again the separation between the stars. 

The short-period binary HD 152248 also presents modulations of the 
X-ray flux that could be reproduced to first order by hydrodynamical 
modelling (see Fig. [2H] and Sana et al. 2004). 

Like these two systems, the binary 7^ Vel (WC8+07.5III) displays 
phase-locked variations of its X-ray fiux. However, these variations are 
simply due to a changing opacity: when the shock cone around the 
0-star is in the line-of-sight, the absorption is much smaller than at 
other phases when the dense wind of the Wolf-Rayet star absorbs most 
of the X-ray fiux (Willis et al. 1995). 
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Figure 27. Left: Schematic view of the HD 93403 bi- 
nary at the time of the four ^MM-Newton observations. 
The top arrows indicate the direction of the observer's 
hne of sight projected on the orbital plane and the dashed 
and the dotted circles correspond to the surfaces of opti- 
cal depth unity for the primary wind at 0.5 keV and 1.0 
keV respectively. Right : X-ray lightcurve of HD 93403 
in the 1.0-2.5 (medium, up), 0.5-1.0 (soft, second from 
above) and 2.5-10.0 keV (hard, third) energy bands with 
1-0" error bars. The last panel yields the relative orbital 
separation between the components of the system while 
the lower panel provides the position angle of the bi- 
nary axis (0° corresponding to the primary star being 
"in front" of the secondary). (From Rauw et al. 2002) 

Figure 28. Left: Schematic view of HD 152248 at the 
time of the six ^MM-Newton pointings. The primary 
star is in dark grey while the secondary is represented in 
light grey. Arrows at left-hand indicate the projection, 
on the orbital plane, of the line-of-sight. Right : Com- 
parison of the observed X-ray luminosities and the results 
from the hydrodynamical simulations: filled squares rep- 
resent the dereddened luminosities of the interaction re- 
gion as predicted by the model; filled circles are the total 
predicted luminosities (including the expected intrinsic 
contribution from the two components of HD 152248); 
open circles show the observed dereddened luminosities. 
(From Sana et al. 2004) 

Finally, the eclipsing binary CPD-41°7742 (09V+B1-1.5V) also presents 
a modulation of the X-ray flux (see Fig.[29]and Sana et al. 2005). How- 
ever, this is a peculiar case of CW since the secondary star has a very 
weak wind. Therefore, it is expected that the wind of the primary di- 
rectly crashes onto the photosphere of the secondary or that it suffers 
radiative braking^, leading to a wind-wind interaction very close to 
the photosphere. Since the interaction region is close to the secondary 
star, its X-ray emission does not reach us when the secondary star is in 
front (only the back of the star is then observed) or when the primary 
occults its companion (Sana et al. 2005). 



In this case, the wind of the primary star is suddenly stopped because of the 
radiative pressure from its companion. The main consequence of this braking is 
to push the contact surface away from the secondary star, preventing the wind to 
actually reach the secondary surface 
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Figure 29. Left: Two schematic views of the geometri- 
cal wind-photosphere interaction model in CPD-41°7742: 
the (a) panel shows a view from above the orbital plane 
whereas the (b) panel presents a view along the line of 
sight just before/after an eclipse. The darker the shad- 
ing of the secondary surface, the higher the X-ray lumi- 
nosity. Right : Predicted unabsorbed flux emitted by 
the wind-photosphere interaction in CPD-41°7742 (thick 
line) compared to the ^MM-Newton observations (filled 
triangles are for MOSl, open squares for M0S2). The 
dashed line gives the intrinsic contribution from the two 
stellar components. (From Sana et al. 2005) 

Figure 30. Aperture masking interferometry has en- 
abled to detect a spiral-shape emitting region around the 
Wolf-Rayet star WR104 (Tuthill et al. 1999). This re- 
gion corresponds to dust forming in the CW shock and 
its shape can be explained as a combination of the or- 
bital rotation and the outward motion of the wind (like 
for a lawn sprinkler), (c) Keck Obs. 

Figure 31. VIA maps of the WR147 system 
(WN8+B0.5V) at 3.6 cm. The asterisks mark the po- 
sition of the two stars. The northern radio emission is 
clearly elongated and not centered on the star: it corre- 
sponds to the radio emission from the CW region. (From 
Contreras & Rodriguez 1999) 

3.4. Other Wavelength Ranges. Colliding winds also have an im- 
pact on the spectrum at other wavelengths. When binaries composed 
of certain Wolf-Rayet stars are observed in the infrared (IR), for ex- 
ample, the formation of dust in the CW region can be detected. In 
eccentric systems, this dust formation is recurrent, since it appears 
only at specific orbital phases, and it leads to a modulation of the in- 
frared emission (e.g. WR140, Williams et al. 1990). Another example 
is WR104 (WC9+lateO-earlyB), a CW binary (CWB) surrounded by 
an IR pinwheel nebula (see Fig. [HUl) that rotates in harmony with the 
orbital period. 

Massive stars generally emit radio waves because of thermal free-free 
transitions in their winds. Such an emission is of the form Fj, oc z/°'^. 
For some binaries, however, an additional, non-thermal emission can be 
observed. This radio emission is of the synchrotron type and is linked 
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to the motion of relativistic electrons in the stellar magnetic field. The 
acceleration of electrons to relativistic velocities is probably achieved 
through the first order Fermi mechanism. Consider a shock in a wind 
moving with a velocity V. A high energy particle crossing the shock 
hardly notices it. But the downstream gas leaves the shock with a lower 
velocity (V/A) then the velocity at which the upstream gas enters the 
shock {V), and the particle gains on average 0.5V/ c by simply crossing 
the shock downwards. The same particle can then be scattered back 
upstream by turbulence without any energy loss. A succession of such 
crossings can thus accelerate the particle to relativistic energies. How- 
ever, the existence of the pre-existing high energy particles, although 
necessary, is not yet explained theoretically. On the other hand, it has 
been proven recently that only shocks between two winds are strong 
enough to accelerate the particles. This implies that all massive non- 
thermal emitters are binaries, and this was also proven recently by 
careful multiwavelength observations (see De Becker, 2005). Note that 
the radio emission arises in the outer regions of the wind, contrary to 
what is seen in X-rays, because the inner parts of the winds are opaque 
to radio waves. The radio emission from CW regions has even been 
imaged in some cases (see Fig. | 



Finally, the relativistic electrons can also give rise to non-thermal X- 
and 7-rays, in addition to synchrotron radio radiation. In this case, 
inverse Compton scattering forces the relativistic electrons to give up 
part of their energy to ambient photons (e.g. the numerous stellar UV 
photons). These UV photons thus become high-energy radiation. In 
addition, relativistic protons can produce neutral pions when they in- 
teract with the ions in the densest regions of the wind: these pions 
subsequently decay into 7-rays. However, this type of emission has not 
yet been detected with certainty. 
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